JOURNAL OF RESEARCH of the National Bureau of Standards— D. Radio Propagation 
Vol. 65D, No. 6, November-December 1961 

Attenuation Coefficients for Propagation at Very Low 
Frequencies (VLF) During a Sudden Ionospheric Dis- 
turbance (SID) ' 

E. T. Pierce 

Contribution from Stanford Research Instituta. Menlo Park, Calif. 

(Received May 11, 1961) 

Attenuation coefficients for propagation under SID conditions are deduced, from records 
of atmospheric noise, for the frequency range of 3.5 to 50 kilocycles per second. These are 
compared with the values under normal daytime circumstances. It is shown that the 
advent of SID implies little change in attenuation between about 12 and 20 kilocycles 
per second; above this range there is a markedly decreased attenuation, while below 12 
kilocycles per second there is an even more pronounced increase in the attenuation coef- 
ficients. The principle of an improved recorder of the changes in atmospheric noise associ- 
ated with SID conditions is outlined; this recorder would discriminate between source 
effects independent of the SID and propagation influences solely attributable to the SID. 



1. Introduction 

The behavior, during propagation under normal 
daytime conditions, of radio waves of very low 
frequency (VLF), is now well established. It is 
known that the propagation is within the wave- 
guide formed by the earth and the lower edge of 
the ionospheric D layer. Attenuation coefficients 
in the VLF range (3 to 30 kc/s) have been deduced 
using the mode theory [Wait, 1957], and have also 
been derived experimentally from observations of 
atmospherics [Chapman and Macario, 1956] and of 
longwave transmitters [Watt and Plush, 1959]. 
The agreement between experimental and theoreti- 
cal results has been good, and the increasing rec- 
gnition of such refinements as the influence of 
round conductivity, the effect of the earth's mag- 
etic field, and so on, is likely to lead to still better 
ccord. 
One influence that considerably affects daytime 
YLF propagation is the sudden ionospheric dis- 
turbance (SID). It is now generally agreed that 
a SID is caused by an increased flux of X-rays 
Emitted by the sun during a solar flare [Friedman, 
1959]. These X-rays enhance the ionization in 
the D layer. One aspect of this increased ionization 
Is that an additional shallow bulge of electrons 
forms upon the underside of the normal D region. 
The exact shape of the bulge depends upon the 
spectral distribution of the X-rays in the incident 
flux [Pierce and Swift, I960]; in all cases, however, 
the depression of the D layer will be greatest at 
the subsolar point, and there will be appreciable 
affects upon the D region within a radius of at least 
3,000 km from the subsolar point. 

As regards propagation at VLF the influence of a 

SID may be represented as a narrowing of the wave- 

uide and, possibly, a change of the effective con- 
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ductivity at the upper boundary of the guide. Con- 
sequences upon propagation may therefore be 
anticipated and these are indeed found. The main 
effects are sudden phase anomalies (SPA) in which a 
phase-change in the sky wave from a manmade 
transmitter is detected, due to the depression of the 
D reflecting layer; sudden anomalies in field strength 
(SFA) which result from the changed phase relation- 
ship between the skywave and the groundwave, 
again from an artificial transmitter; and sudden 
enhancements of atmospherics (SEA) in which an 
increase of the noise due to distant thunderstorms is 
observed (typically at a frequency of 27 kc/s). SPA 
and SFA studies indicate that for a specific electron 
density an average SID lowers the D layer by some 
5 km [Bracewell and Straker, 1949]. 

Although, as already mentioned, VLF attenuation 
coefficients have been widely obtained for normal 
conditions, the same information does not appear to 
have been derived under SID circumstances. The 
essential raw material, namely the changes observed 
in field strengths, has however been available for 
some time. It is the purpose of this note to attempt 
to obtain some indication of the VLF attenuation 
coefficients, during a SID, from data upon atmos- 
pherics. 

2. Normal Daytime Attenuation Coefficients 

Figure 1 depicts the relationship between the 
attenuation coefficient and frequency during normal 
daytime propagation. Two curves are shown . The 
first curve was deduced by Chapman and Macario 
[1956] from observations of atmospherics. The 
second is that derived by Watt and Plush [1959], 
largely on the basis of observations of manmade 
transmitters, and is described as being an average 
for a mixed path. There are divergencies between 
the two curves but these are not entirely unexpected. 
Even if D layer ionospheric conditions were uniform 
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Figure 1. Attenuation coefficients for normal daytime condi- 
tions. 

during daylight irrespective of solar altitude, such 
influences as ground conductivity and the direction 
of propagation with respect to the earth's magnetic 
field would still be important. These influences 
may lead to systematic deviations from "average" 
behavior. For instance, it is not unusual for the 
main sources of atmospherics received at a particular 
station to lie within a restricted geographical area; 
this implies that the attenuation coefficients, ob- 
tained from such observations of atmospherics, will 
be representative of a certain path, of its orientation 
with regard to the magnetic field and of the con- 
ductivity along its ground boundary, rather than of 
more general conditions. 

3. Basic Data and Its Analysis 

Wait [1958] has given the following expression for 
the vertical electric field E (in millivolts/meter) at a 
great circle distance d from a source radiating P kw: 
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where a is the radius of the earth, h the height of 
the ionospheric reflecting layer in kilometers, X is 
the wavelength in kilometers, and a the attenuation 
factor. Equation (1) assumes that only a single 
mode is dominant, but this is a fair approximation 
up to 30 kc/s and for distances in the range 1,000 
to 2,000 km, and almost completely true for d greater 
than 2,000 km; thus (1) is valid in the study of 
distant atmospherics. Equation (1), however, is 
likely to be increasingly incorrect as the frequency 
rises above 30 kc/s. Under normal conditions 
(suffix N) (1) may be written as 
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If atmospheric noise is being received, then the 
advent of SID conditions (suffix S) will not cause 
P or d to change and 
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Hence 



h N E N /h s E s = exp (a s — a^d, 



(4) 



it being assumed that E N and E s are compared at 
the same frequency. Equation (4) may be con- 
verted into the form 



(a s -a N )d=2XlO* hg l0 (h N E N /h s E s ), 



(5) 



in which the units for a s and a N are now the familiar 
db/1,000 km. Thus if records of E N and E s are 
available together with estimates for d, h Nj and 
hs, (as—ax) can be deduced, and, knowing a N , a s 
obtained. 

Four sets of data in the published literature give 
the ratio E N /E S at different frequencies and are 
therefore particularly suitable for the application of 
eq (5) . These instances are in the papers by Gardner 
[1950], by Obayashi, Fujii, and Kidokoro [1959], 
and by Obayashi [I960]; this last paper contains two 
examples. Internal evidence within Gardner's paper 
suggests that 1,500 km is a good estimate for the 
distance of the noise sources; the reliability of this 
estimate is enhanced by the fact that Gardner 
worked in close collaboration with the Sferics 
Locating Network of the British Meteorological 
Office. It remains to determine h N /h s before (a s — 
a N ) can be calculated. h N jh s is not far removed 
from unity and the errors in its estimation are 
considerably less than those implicit in assigning 
a value to d. 70 km is a commonly accepted value 
for h N , while h s may be taken as (65 ±5) km [Brace- 
well and Straker, 1949] the extremes corresponding 
to the slightest and to the most intense SID. Thus 
W^s— 1.08. Values of E N /E S were obtained from 
Gardner's paper and, using eq (5), (a 5 -%)was 
calculated over the frequency range 3.5 to 50 kc/s. 
The results are represented by curve (a) on figure 2. 

Similar information to that of Gardner is given 
in the much later paper by Obayashi and coworkers 
[1959]. Here the frequency range is restricted to 
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Figure 2. Difference between attenuation coefficients for 
normal daytime and for SID conditions. 
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10 to 50 kc/s. There is no direct evidence in the 
paper as to the distances of the sources of atmos- 
pheric noise. However, the SID concerned occurred 
between 0430 and 0500 U.T. on 16 August 1958. 
Atmospheric noise records for this day indicate no 
local sources, while sferics hearing observations 
covering the period of the SID show the activity 
to be to the east and to the southwest of Japan 
[IGY Data, Japan, 1959]. It seems plausible to 
associate these directions respectively with the 
main storm centers in the Americas and in SE Asia. 
The activity in the Americas will be dying with the 
advent of local night; that in SE Asia will be rapidly 
increasing the time being adjacent to local noon. 
In atmospheric noise measurements the nearest 
storms are the dominant factors [Horner, I960]; 
thus the SE Asia center is likely to be more important 
than that in the Americas. Furthermore, it is the 
closer parts of the SE Asia center that will be most 
vital; this effect is enhanced by the influence of time 
since maximum thundery activity occurs in the local 
afternoon, and it is generally in the nearer parts of 
the SE Asia center to Japan, that, at the time of 
the SID, the local maximum is most closely ap- 
proached. Considering all factors, it seems reason- 
able to identify the principal sources of the atmos- 
pherics as being within the nearest landmasses to 
the southeast of Japan, and to take d as being 2,000 
km. Using h N /h s =l.0S, the values of (0^—0^) 
were calculated as before and are represented by 
curve (b) on figure 2. 

Obayashi [1960] has given two further examples 
of enhancements of atmospheric noise which are 
suitable for analysis. These are for 29 March 1959 
and 16 June 1959. The distance of the sources 
of noise for the first instance is estimated at 2,000 
km, and Obayashi states that the storms were 
rather closer for the second enhancement; d= 1,500 
km would seem to be a reasonable mean estimate in 
this case. Using h N /h s =l.08 as before, the results 
obtained are shown in curves (c) and (d) on figure 2. 

The agreement between the four curves in figure 
2, for the range of 8 to 35 kc/s, is surprisingly good 
considering the speculative nature of the estimates 
for d in particular. Above 35 kc/s and below 8 
kc/s, the signal strengths involved are low, and 
thus the results are more prone to errors. This 
effect may account for the spread between 35 and 
50 kc/s, but seems quite inadequate to explain the 
discrepancies below 8 kc/s. Here the divergence 
between Gardner's and Obayashi's work is very 
considerable. The estimate for d is probably better 
for Gardner's observations than for Obayashi's; 
however, this factor is insufficient to enable the 
disagreements to be reconciled, and further experi- 
mental work is obviously desirable. 

An intriguing feature of three of the curves on 
figure 2 is a maximum value and then a drop in 
{a s — a N ) as frequency decreases below 10 kc/s. If 
this trend is continued it is probable that at some 
lower frequency, possibly in the extremely low fre- 
quency (ELF) range «3 kc/s), (a s — a N ) will become 
negative as it is at frequencies above about 15 kc/s. 



It might be predicted that a SID will produce an 
enhancement of atmospheric noise above 15 kc/s; a 
decrease at intermediate frequencies; and an en- 
hancement again in the ELF band. This kind of 
predicted behavior, which would have analogies with 
the usual variation of a N with frequency, could 
readily be verified by an extension of experimental 
observations to lower frequencies. 

4. Discussion 

Figure 3 is an attempt to summarize the preceding 
information. Curve (a) is an average representation 
for a N derived from the two graphs of figure 1. 
Curve (b) is a composite graph for a s ; it was ob- 
tained by considering the four curves of figure 2, 
deducing weighted mean values for (a N — a s ), and 
using these, in conjunction with curve (a) of figure 3 
to derive curve (b). 

It is apparent from figure 3 that the frequency 
of 27 kc/s which has been established for many 
years [Bureau, 1937] as an indicator of enhance- 
ments in atmospherics is a good choice, since the 
minimum of a s is close to 27 kc/s and the change 
between <x N and a s is considerable. It is equally 
true that the differences between a N and a s are 
most marked for the lower frequencies, and the 
decrease in atmospherics over that frequency range 
is potentially a more sensitive indicator of the effects 
of a SID than is the enhancement at frequencies 
exceeding some 15 kc/s. This is particularly so 
if the normal storm centers are sufficiently close 
to the observing station for adequate daytime 
signal strength to be always present. An interest- 
ing feature of figure 3 is that the band (15 to 20 
kc/s) occupied by the most important VLF trans- 
mitters is little affected in attenuation by the onset 
of SID conditions. 

The normal atmospherics receiver on 27 kc/s 
which is usually intended as an indicator of the 
sudden enhancements associated with a SID, can 
also give increased signals produced by source 
effects and not by the propagation changes due 
to the SID. Three main causes may be assigned 
to such increases, namely, an augmentation of the 



(a) NORMAL DAYTIME 

-- (b) SID. 
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Figure 3. Composite values for attenuation coefficients under 
normal daytime and under SID conditions. 
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activity in an existing storm center; the develop- 
ment of a new storm center; and the movement 
of existing activity closer to the observing station. 
It seems possible, upon the basis of figure 3, to pro- 
duce a reciprocal instrument which would not 
respond to the above source effects, but would 
only register a change in propagation conditions 
such as that produced by a SID. The principle 
would be to select two frequencies / x and jf 2 , one on 
each side of the normal minimum for a N at around 
18 kc/s, and for which the values of a N are equal; 
as an example figure 3 indicates that 9 and 35 kc/s 
could be chosen for j x and j 2 . Then if E(f 1 ) and 
E(J 2 ) are the measured fields at j x and / 2 , the ratio 
E{ji)jE{j 2 ) should remain effectively constant irre- 
spective of alterations in the strengths of the sources 
of atmospheric noise, or of changes in the distances 
of the sources, provided that these are not such as 
greatly to modify the validity of the basic one- 
mode assumption; this implies distances always 
greater than some 1,000 km. However if prop- 
agation conditions are affected as in a SID, the 
ratio E(f 1 )/E(f 2 ) will alter very considerably be- 
cause of the changes in attenuation coefficient 
from a N to a s . The effect will be especially marked 
being reciprocal at the two frequencies, as is shown 
in figure 3 for the selected frequencies of 9 and 35 
kc/s. Construction of a simple recorder presenting 
the ratio E(f 1 )/E(f 2 ) directly offers no great difficulty. 
It is outside the scope of this paper to examine the 
interpretation of the a s values of figure 3 upon the 
waveguide mode theory. However, there is no 
reason to believe that the values are incompatible 
with the theory, indeed the opposite is rather the 
case and the results are further evidence — although 
this is hardly necessary — of the general validity of 
the mode interpretation. For instance, it can 
readily be deduced from the work of Wait [1957] 
that a reduction in h from 70 to 60 km, without any 



considerable associated change in ionospheric con- 
ductivity, produces two graphs very similar to those 
depicted on figure 3. In each case, there is a recipro- 
cal relative displacement of the two curves, pivoting 
around a common point. 
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